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Accomplishments: 

A.  The  Relationship  Between  Fock-Space  Multi- reference  CC  and  EOM-CC: 

This  year  much  effort  was  devoted  to  detailing  the  relationships  between  two 
different  methods  of  treating  excited  states  in  CC  theory  [9)12, 14, 17,18].  Today  the 
EOM-CC  approach  is  well-established,  after  being  introduced  by  us  in  the  90's.  It  is 
widely  used  for  many  problems  like  the  treatment  of  electronic  (UV-vis)  and 
photoelectronic  spectra.  Even  more  significantly,  it  is  quite  valuable  as  a  ‘target  state 
method',  meaning  the  state  of  interest  can  be  better  described  via  the  EOM-CC 
framework  than  it  can  be  within  the  usual  single-reference  CC  approach.  An 
example  is  the  treatment  of  open-shell  singlet  states  that  are  formally  exact  in  the 
operationally  single-reference  CC  approach  (ie  no  decisions  except  choice  of  basis 
and  level  of  correlation),  but  would  otherwise  require  the  use  of  two  determinants  (or 
a  multi-reference  description)  in  a  more  conventional  approach.  This  occurs  because 
EOM-CC  diagonalizes  a  matrix  that  would  have  both,  equivalent  determinants  in  it, 
allowing  them  to  be  weighted  however  warranted.  The  same  is  true  for  more 
complicated  low-spin  cases.  Even  more  complicated  multi-reference  character  can 
be  introduced  by  using  the  flexibility  inherent  in  EOM-CC  to  start  from  a  problem  with 
a  different  number  of  electrons  to  correctly  describe  a  problem,  as  in  using  O3'2  to 
obtain  a  description  of  O3,  by  removing  two  electrons  from  the  anion  to  correctly 
describe  the  quasi-degeneracy  between  the  two  homo  orbitals  in  ozone. 

The  alternative  FS  approach  uses  a  multi-reference  description  to  accomplish  some 
of  the  same  things,  but  has  the  advantage  over  EOM-CCSD,  eg,  in  that  that  method 
would  require  the  diagonalization  of  a  matrix  of  the  dimension  of  all  single  and 
doubleexcitations,  while  the  FS  approach  would  get  its  answers  from  only  the 
diagonalization  of  a  matrix  of  the  dimension  of  singles  alone,  which  are  much  fewer. 
Furthermore,  the  FS  approach  formally  derives  from  an  exponential  ansatz  for  the 
excited  states,  while  the  EOM-CC  excited  states,  contrary  to  the  ground  state,  use  a 
linear  (Cl-like)  operator.  The  numerical  distinction  is  modest,  but  can  be  observed  in 
the  limit  where  an  excited  state  would  separate  into  A+  and  B-.  We  refer  to  this  as 
‘charge-transfer’  separability.  All  other  excited  states  would  go  to  A*  +  B*  in  either 
method. 

Hence,  in  an  attempt  to  get  the  best  of  both  worlds,  we  have  embarked  upon  the 
detailed  relationship  between  the  two  methods.  This  has  been  accomplished  by 


invoking  the  Intermediate  Hamiltonian  approach  in  FS-CC,  and  relating  it  through  the 
Bloch  equation  to  EOM-CC.  Because  of  the  formal  relationships  we  have 
established,  we  now  have  a  route  toward  treating  any  sector  of  FS,  ie  any  number  of 
electrons  outside  a  closed  shell,  which  had  previously  proved  to  be  almost 
impossible  beyond  the  1  hole  1  particle  sector;  by  straight-forward  matrix 
diagonalization,  a  very  stable  procedure.  Also,  EOM-CC  has  always  provided  a 
route  toward  treating  properties,  but  FS-CC  has  not.  Exploiting  the  inter¬ 
relationships,  we  now  can  overcome  this  limitation  of  the  FS-CC  approach  and  gain 
from  its  simpler  computational  aspect  and  superior  formal  aspects. 

In  the  course  of  doing  this,  we  also  presented  the  first  results  for  FS-CCSDT[14], 
and,  in  particular,  for  the  FS-CCSDT  for  triplet  states[17,18].  We  also  investigated 
the  route  to  low-spin  states  in  EOM-CC  via  the  spin-flip  idea  of  Krylov,  meaning  use 
the  high-spin  triplet  or  higher  as  a  reference  and  flip  the  unpaired  spins  to  get  an 
approximation  to  a  low-spin  eigenstate  via  the  EOM-CC  approach[12,  23].  This 
further  enhances  what  states  can  be  described  within  EOM-CC  with  minimal 
decisions  for  the  user. 

B.  A  New  Hierarchy  for  CC  theory: 

In  last  year’s  report,  we  focused  on  this  development.  Briefly,  in  the  nCC 
Heirarchy[2],  we  insist  that  the  method  be  exact  for  2,  nCCSD,  3,  nCCSDT,  and  4, 
nCCSDTQ  electrons.  Hence, 

E(CISD)=E(CCSD)=E(nCCSD)  for  n=2; 

E(CISDT)=E(CCSDT)=E(nCCSDT)  for  3, 
E(CISDTQ)=E(CCSDTQ)=E{nCCSDTQ)  for  4. 

All  nCC  methods  are  also  extensive.  Hence,  whereas  when  we  obtain  the  energy 
from  Cl  we  have  the  well-known  size-extensivity  failure, 

E[CISD(A+B)]^ECISD(A)+ECISD{B)  when  A  and  B  are  non-interacting  two  electron 
units,  it  follows  that 

E[CISDTQ]  =  E[CCSDTQ]=E[nCCSDTQ] 

=  E[CCSD(A)]+E[CCSD(B)] 

=  E[nCCSD{A)]+E[nCCSD(B)] 

That  is,  all  products  of  non-interacting  two-electron  units  are  exact  in  nCCSD. 
Similarly,  all  products  of  2,  3  and  4  electrons  also  give  the  exact  energy  for  all  their 
combinations.  Hence,  nCC  enables  us  to  think  about  electron  correlation  in 
conceptually  appealing  2,  3,  and  4  electron  units.  This  also  facilitates  our 
understanding  of  bond  breaking,  as  such  units  will  occur  as  the  separated  products 
after  a  bond  fission.  If  we  can  ensure  a  smooth  transition  to  the  products  by  CC, 
then  we  will  have  a  route  toward  improving  the  bond  breaking  in  any  single 
reference  (UHF,  RHF,  or  other  like  Brueckner  or  a  KSDFT  reference. 


Because  of  nCCSD  being  exact  for  two  eiectrons,  and  all  products  of  two  electrons, 
it  is  the  natural  zeroth-order  approximation  for  the  separated  electron  pair  bonds  in 
chemistry,  which  we  are  further  studying  in  our  Natrual  Linear  Scaling  CC  approach 
(NLSCC)[10,11] .  Furthermore,  the  transition  from  CCSD  to  nCCSD  eliminates  the 
most  time-consuming  steps  in  the  evaluation  of  the  non-linear  terms  in  CCSD.  The 
normal  scaling  of  such  terms  is  ~N3  with  the  number  virtual  orbitals.  However,  they 
can  be  eliminated  leaving  just  the  non-linear  terms  that  scale  as  ~N2.Nonetheless, 
the  method  is  exact  for  two  electrons,  and  all  products  of  two  electrons  as  it  is  still 
composed  solely  of  linked  diagrams  as  required  by  its  extensivity.  Furthermore,  it 
still  retains  all  the  invariance  properties  of  CCSD,  meaning  that  it  is  invariant  to 
transformations  among  just  occupied  or  virtual  orbitals.  This  distinguishes  this 
method  from  CEPA-like  approximations  that  are  not  invariant  to  orbital  rotations. 
Finally,  nCCSD  has  the  same  orbital  insensitivity  as  does  CCSD,  since  the  exp(Ti) 
operator  is  not  truncated.  The  nCC  approach  coupled-with  NLSCC  is  being 
implemented  into  ACES  III,  the  completely  newly  written,  massively  parallel  version 
of  ACES  II.  In  this  way  very  large  molecules  will  be  studied. 

C.  ACES  III  and  massively  parallel  implementations. 

From  scratch,  we  have  written  ACES  III  over  the  last  few  years  [15,  20]  and  are  now 
starting  to  reap  the  benefits.  As  shown  below,  we  have  achieved  virtually  perfect 
scaling  to  1024  processors  and  even  85%  scaling  up  to  4096.  We  have  also 
successfully  run  on  >30000  processors.  What  this  means  is  that  soon  it  will  be 
almost  routine  to  make  CC  applications  on  10,000  processors,  which  is  normally 
considered  the  terascale,  with  the  petascale  in  prospect.  Furthermore,  since  ACES 
III  provides  CCSD  and  CCSD(T),  for  open  and  closed  shells,  with  any  single 
determinant  reference  function,  and  analytical  gradients  for  same,  [plus  EOM-CC  for 
excited  states,  [lus  its  IP-  and  EA-EOM  variants,  it  can  be  used  in  plethora  of 
applications  to  provide  ‘benchmark’  accuracy  that  were  previously  thought  to  be  well 
beyond  the  reach  of  CC  theory.  This  development  will  offer  a  paradigm  change  in 
the  field. 


CCSO  scaling  results  for  various  molecules 


D.  Further  progress  in  the  development  of  ab  initio  dft 

The  idea  of  ab  initio  dft  is  that  the  Kohn-Sham  density  functional  one-electron 
structure  can  be  built  upon  ab  initio  correlated  coupled -cluster  and  many- body 
perturbation  theory  methods,  to  correct  the  principal  failing  of  dft  methods;  namely 
that  they  do  not  converge  to  the  right  answer  in  the  basis  set  and  correlation  limit, 
unlike  ab  initio  wavefunction  methods.  On  the  other  hand,  wavef unction  methods 
require  a  two-particie  structure,  which  limits  their  application  to  large  molecules, 
unlike  DFT.  Progress  can  be  made  by  combining  the  best  of  both  worlds.  We  have 
previously  shown  a  number  of  results  that  conclusively  prove  the  formal  and 
numerical  superiority  of  such  an  approach,  but  always  at  a  cost  in  its  implementation 
compared  to  conventional  DFT.  Papers  [4,5,13]  further  contribute  to  this 
development. 


Paper  [5]  addresses  the  issue  of  quasi-degeneracy  in  ab  initio  dft,  using  the  Be 
isoelectronic  sequence  as  an  example.  For  this  sequence,  previous  attempts  by 
other  authors  to  obtain  the  correct  exchange  and  correlation  potentials  suffer  from 
divergence.  For  our  ab  initio  dft  approach,  this  does  not  happen  because  of  the 
superior  choice  we  make  for  our  unperturbed  hamiltonian,  compared  to  the  usual 
Kohn-Sham  choice  that  is  invariably  used  by  DFT  practioners.  However,  for  the  Be 
sequence,  there  is  still  a  need  to  incorporate  higher-order  effects  than  is  possible  by 
using  a  MBPT(2)  functional.  Hence,  we  introduced  a  new  coupled-cluster  based 


perturbation  theory  to  add  higher-order  effects  into  what  is  still  an  attractive  second- 
order  approach.  We  call  the  new  approach,  CCPT2.  With  it  we  are  able  to  show 
much  better  converged  results  for  the  Be  sequence.  This  approach  also  has  some 
advantages  for  other  coupled-cluster  problems,  which  we  are  pursuing. 

Paper  [13]  asks  the  interesting  question  can  we  use  a  fraction  of  HF  exchange  to 
make  the  correlation  potential  that  ab  initio  dft  creates  to  be  obtained  more  easily, 
and  can  we  define  an  exact  one-particle  theory  such  that  the  orbital  energies 
obtained  from  the  one-particle  Hamiltonian  will  provide  virtually  exact  ionization 
potentials  for  all  the  occupied  states  of  a  molecule.  We  can  indeed  do  that  with  the 
simple  expedient  of  taking  50%  HF  exchange  together  with  50%  OEP  exchange,  but 
it  is  necessary  to  also  use  our  usual  semi-canonical  OEP-MBPT2  correlation 
potential  to  make  this  happen.  This  is  a  fascinating  result  as  it  means  that 
numerically  one  can  invent  a  one-particle,  Hamiltonian  that  gives  exact  results.  So 
this  can  be  viewed  as  a  frequency  independent  approximation  to  the  Dyson  operator 
or  as  a  generalization  of  DFT  that  insists  upon  this  IP  property  as  an  exact  condition 
instead  of  the  density  that  defines  KS  theory. 

E.  Our  intent  in  new  work  is  to  introduce  an  alternative,  exact,  correlated  one- 
particle  theory  for  electronic  structure,  for  both  its  computational  aspects  and  its 
conceptual  ones.  It  will  correctly  describe  self-interaction  and  charge-transfer,  and 
dispersion.  It  also  provides  a  systematic  route  toward  exact  solutions  as  more 
particles  are  added.  We  will  also  address  the  question  of  what  quantities  can  we 
expect  to  be  able  to  describe  correctly  within  such  a  one-particle,  correlated  orbital 
theory  (COT).That  is,  we  insist  that  an  hef,{1 )  be  created  such  that 

hef,(1)tPi(1)=Ljtpi(1),  l,=-l(i) 
he,f(1)<pa(1)=na<pa(1),  □a=A(a) 


In  each  case,  the  eigenvalues  are  the  exact,  principal  ionization  potentials,  l(i),  and 
electron  affinities,  A(a).  The  meaning  of  'principal'  means  the  energies  required  to 
remove  an  electron  from  any  of  the  iA{th}  occupied  orbitals,  or  to  add  an  electron  into 
any  of  the  aA{th}  unoccupied  orbitals.  For  the  latter,  many  such  energies  will  be 
negative,  that  is,  will  not  correspond  to  bound  anions  and  observed  electron  affinities. 
For  I's,  most  of  the  principal  ones  will  usually  be  seen  experimentally  except  when  a 
'shake-up'  has  a  very  large  intensity  that  masks  the  principal  in  that  energy  range.  For 
our  purposes,  a  shake-up  in  not  a  one-particle  effect,  so  it  falls  outside  the  scope  of 
our  present  exact  COT.  If  we  do  an  IP-EOM-CC  calculation,  matrix  diagonalization 
will  provide  a  set  of  n  principal  eigenvalues  and  vectors  for  n-electrons  (orbitals). 
Similarly,  EA-EOM-CC  will  provide  the  same  A's  for  m-n  unoccupied  orbitals,  where 
m  is  the  dimension  of  the  MO  basis  set.  These,  of  course,  will  be  exact  in  the  limit  of 
electron  correlation  and  basis  set.  They  are  also  equivalent  to  those  obtained  from 
effective  one-particle  theories  as  we  can  prove.  Because  he,f  is  non-Hermitian,  the 


orbitals  will  compose  both  a  right-hand  and  a  left-hand  set,  with  the  same 
eigenvalues. 

Hence,  we  first  establish  that  an  'exact'  one-particle  theory  of  the  type  we  propose 
exists.  We  achieve  this  by  construction [25].  Furthermore,  those  equations  can  always 
be  evaluated  by  doing  coupled-cluster  calculations.  This  is  notably  contrary  to  DFT, 
which  is  not  based  on  'constructive'  proofs,  but  instead  on  existence  proofs.  Hence, 
there  is  seldom  any  independent  way  to  assess  the  veracity  of  chosen 
approximations.  But  in  our  new  theory  an  exciting  prospect  is  that  the  equations  can 
be  used  independently  of  doing  the  actual  CC  calculations,  provided  that  an 
adequate  approximation  for  the  two-partcle  interactions  can  be  made,  like  that  due  to 
Colie  and  Salvetti,  or  for  the  two-particle  density  matrix  in  general  as  in  density  matrix 
functional  theory.  Unlike  DFT,  such  an  approach  provides  a  litmus  test  for  any  two- 
particle  approximation  considered,  since  the  eigenvalues  of  the  associated  potential 
have  to  reflect  the  exact  properties  we  establish.  Furthermore,  in  our  constructive 
method,  we  can  always  make  relevant  coupled-cluster  calculations  to  see  what  the 
right  answers  ought  to  be,  and  use  this  to  immediately  assess  various 
approximations.  We  are  currently  making  the  initial  computations  to  illustrate  its  use. 
Among  other  things,  this  theory  will  introduce  a  set  of  ionization(l)  and  electron 
attached  (A)  orbitals  that  associate  with  each  orbital  the  exact  principal  I  and  A 
energies.  So  the  difference  between  the  HOMO  and  LUMO  levels  will  be  exactly  the 
smallest  l-A.  In  the  case  of  solids  this  is  the  band  gap,  while  for  molecules,  this  theory 
should  offer  a  new  conceptual  (and  computational)  basis  for  MO  theory  based  upon 
exact  I  and  A  values. 

Furthermore,  I  would  also  suggest  that  having  this  energy  based  scheme  rather  than 
a  density  based  one  (DFT)  is  vastly  more  suited  to  photo-electron  spectra  in  general, 
a  pervasive  component  of  AFOSR’  programs,  such  as  in  its  applications  to  ionic 
liquids,  and  in  the  identification  and  classification  of  clusters  via  such  PES 
experiments  using  anions,  cations,  and  neutrals.  The  new  theory  I  am  proposing  has 
to  get  I's  and  A's  right,  unlike  DFT.  And  when  needed,  the  corresponding  CC 
applications  can  be  made  to  keep  the  simplified  theory,  which  is  geared  toward  very 
large  molecule  applications,  on  target. 
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January  2007  -  National  Chemical  Laboratory,  Pune,  India 


Changes  in  research  objectives,  if  any:  None 
Change  in  AFOSR  program  manager,  if  any:  None 
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